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Abstract 



Correlated 5 distributions can be used to describe the clutter seen in images obtained with coherent illumination, 
as is the case of B-scan ultrasound, laser, sonar and synthetic aperture radar (SAR) imagery. These distributions 
are derived using the square root of the generalized inverse Gaussian distribution for the amplitude backscatter 
within the multiplicative model. A two-parameters particular case of the amplitude 9 distribution, called 
3^4, constitutes a modeling improvement with respect to the widespread Xa distribution when fitting urban, 
forested and deforested areas in remote sensing data. This article deals with the modeling and the simulation 
of correlated S^-distributed random fields. It is accomplished by means of the Inverse Transform method, 
applied to Gaussian random fields with spatial correlation. The main feature of this approach is its generality, 
since it allows the introduction of negative correlation values in the resulting process, necessary for the proper 
explanation of the shadowing effect in many SAR images. 
Keywords: image modeling, simulation, spatial correlation, speckle. 



1 Introduction 



The demand for exhaustive and controlled clutter measurements in all scenarios would be alleviated if plausible 
data could be obtained by computer simulation. Clutter simulation is an important element in the development 
of target detection algorithms for radar, sonar, ultrasound and laser imaging systems. Using simulated data, 
the accuracy of clutter models may be assessed and the performance of target detection algorithms may be 
quantified with controlled clutter backgrounds. This article is concerned with the simulation of random clutter 
having appropriate both first and second order statistical properties. 

The use of correlation in clutter models is significant and relevant since the correlation effects within the 
clutter often dominate system performance. Models merely based on single-point statistics could, therefore, 
produce misleading results, and several commonly used forms for clutter statistics fall into this category. 

The statistical properties of heterogeneous clutter returned by Synthetic Aperture Radar (SAR) sensors 
have been largely investigated in the literature. A theoretical model widely adopted for these images assumes 
that the value in every pixel is the observation of an uncorrelated stochastic process Za, characterized by single- 
point (first order) statistics. A general agreement has been reached that amplitude fields are well explained 
by the 3Ca distribution. Such distribution arises when coherent radiation is scattered by a surface having 
Gamma-distributed cross-section fluctuations. Though agricultural fields and woodland are very well fitted by 
this distribution, it is also known that it fails giving accurate statistical description of extremely heterogeneous 
data, such as urban areas and forest growing on undulated relief. 

As discussed in [UI2], another distribution, the Sa law, can be used to describe those extremely heterogeneous 
regions, with the advantage that it has the 3Ca distribution as a particular case. This distribution arises in all 
coherent imaging applications as a result of the action of multiplicative speckle noise on an underlying square 
root of a generalized inverse Gaussian distribution. The main drawback of this general model is that it requires 
an extra parameter, besides its theoretical complexity. 

Nevertheless, it can be seen in HJ \5\ that a special case of the Sa distribution, namely the Sa l aw > 
which has as many parameters as the %a distribution, is able to model with accuracy every type of clutter. 
As a consequence, efforts have been directed toward the simulation of Sa textures, but no exact method for 
generating patterns with arbitrary spatial autocorrelation functions has been envisaged so far, in spite of it 
being more tractable than the %a distribution. 

As previously stated, spatial correlation is needed in order to increase the adequacy of the model to real 
situations. This paper tackles the problem of simulating correlated 5 a fields. 



The main properties and definitions of the S\ clutter are presented in this section, starting with the first order 
properties of the distribution and concluding with the definition of a Sa stochastic process that will describe 
Za fields. 

2.1 Marginal properties 

The S^(a,7,n) distribution is characterized by the following probability density function: 



being n > 1 the number of looks of the image, which is controlled at the image generation process, and It(0 
the indicator function of the set T. The parameter a describes the roughness, being small values (say a < —15) 
usually associated to homogeneous targets, like pasture, values ranging in the (—15, —5] interval usually observed 



2 Correlated $ A clutter 



fz A (z, (a,7,n)) 




•!(o,+oo)(2), a<0,7>0, 



(1) 
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in heterogeneous clutter, like forests, and big values (—5 < a < for instance) commonly seen when extremely 
heterogeneous areas are imaged. The parameter 7 is related to the scale, in the sense that if Z is S5i( a ; E n ) 
distributed then Za = y^fZ obeys a SV a ' 7' n ) l aw - 

A SAR image over a suburban area of Miinchen, Germany, is shown in Figure [TJ It was obtained with 
E-SAR, an experimental polarimetric airborne sensor operated by the German Aerospace Agency (Deutsches 
Zentrum fur Luft- und Raumfahrt - DLR e. V.) The data here shown were generated in single look format, 
and exhibit the three discussed types of roughness: homogeneous (the dark areas to the middle of the image), 
heterogeneous (the clear area to the left) and extremely heterogeneous (the clear area to the right). 




Figure 1: E-SAR image showing three types of texture. 
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The r-th moments of the 9^(ct,7,n) distribution are 

\nJ I (— a)l (n) 

when — r/2 < a < the r-th order moment is infinite. Using equation @ the mean and variance of a S^(a, 7, n) 
distributed random variable can be computed: 

^r(n+|)r(-a-i) 

V>Za = ' = — 





n r(n)r(— a) 
7 [nT 2 (n){-a - l)T 2 (-a - 1) - T 2 (n + ±)r 2 (-a - ±)] 
^ ~ nr 2 (n)r 2 (-a) 



Figure [5] shows three densities of the S^(a, 7, n) distribution for the single look (n = 1) case. These densities 
are normalized so that the expected value is 1 for every value of the roughness parameter. This is obtained using 
equation ^ for setting the scale parameter 7 = 7 Q ,„ = n (T(— a)T(n)/ (T(— a — l/2)T(n + 1/2))) 2 . These densi- 
ties illustrate the three typical situations described above: homogeneous areas (a — —15, dashes), heterogeneous 
clutter (a = —5, dots) and an extremely heterogeneous target (a = —1.5, solid line). 




Figure 2: Densities of the 5 A (a,j a ,3,3) distribution. 



Following Barndorff-Nielsen and Blaesild [BJ, it is interesting to see these densities as log probability functions, 
particularly because the S A i s closely related to the class of Hyperbolic distributions [7] . Figure [3] shows the 
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densities of the 9^(— 3, 1,1) and N(3tt/16, 1/2 — 97r 2 /256) distributions in semilogarithmic scale, along with 
their mean value /j, = 37r/16. The parameters were chosen so that these distributions have equal mean and 
variance. The different decays of their tails is evident: the former behaves logarithmically, while the latter 
decays quadratically. This behavior ensures the ability of the S*a distribution to model data with extreme 
variability. 




0.5 1.0 1.5 2.0 2.5 5.0 

normalized gray scale 

Figure 3: Densities of the S4 and Gaussian distributions with same mean values [i = 37r/16 in semilogarithmic 
scale. 

Besides being essential for the simulation technique here proposed, cumulative distribution functions are 
needed for carrying out goodness of fit tests and for the proposal of estimators based on order statistics. It 
can be seen in [31 [8, 9J that the cumulative distribution function of a ^ A {a,^j,n) distributed random variable is 
given, for every z > 0, by G(z, (a, 7, n)) = T2 n ,-2a{—cz 2 /j), where T Sjt is the cumulative distribution function 
of a Snedecor's F Si t distributed random variable with s and t degrees of freedom. Both T v and T7. 1 are readily 
available in most platforms for computational statistics. 

The single look case is of particular interest since it describes the noisiest images and it exhibits nice 
analytical properties. The distribution is characterized by the density f(z; a, 7, 1) = — ^§ 2(7 + z 2 ) a_1 I(o, o)(- z ), 
whith —a, 7 > 0. Its cumulative distribution function is given by F(t) = 1 — (l + < 2 /7) Q I(o,oo) (*)> an d its 
inverse, useful for the generation of random deviates and the computation of quantiles, is given by F (t) = 

( 7 ((i-t) 1 /"-i)) 1/2 v 1) W. 
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2.2 Correlated clutter 

Instead of denning the model over Z 2 , in this section a realistic description of finite-sized fields is made. Let 
Za = {ZA(k,£))o<k<N—i.o<e<N-i be the stochastic model that describes the return amplitude image. 

Definition 1 We say that Za is a (a, 7, n) stochastic process with correlation function pz A (in symbols Za 
~ (S5i( a J li n )i Pz A )) if f or a M < i, j, k, £ < N — 1 holds that 

1. ZA(k,£) obeys a S^(a,7,rt) law; 

2. the mean field is \xz A = E(ZA(k,£)); 

3. the variance field is o~ 2 Za = Var{ZA(k, £)); 

4- the correlation function is /Oz A ((i,j), (k,£)) — (E(Za(i, j)ZA(k,£)) — Hz A ) l a z A - 

The scale property of the parameter 7 implies that correlation function pz A and 7 are unrelated and, 
therefore, it is enough to generate a Z\ ~ (S^(oi, l,n), pz A ) field and then simply multiply every outcome by 
7 1 / 2 to get the desired field. 

This paper presents a variation of a method used for simulation of correlated Gamma variables, called 
Transformation Method, that can be found in [10] . This method can be summarized in the following three 
steps: 

1. Generate independent outcomes from a convenient distribution. 

2. Introduce correlation in these data. 

3. Transform the correlated observations into data with the desired marginal properties |llj . 

The transformation that guarantees the validity of this procedure is obtained from the cumulative distribu- 
tion functions of the data obtained in step [21 and from the desired set of distributions. 

Recall that if U is a continuous random variable with cumulative distribution function Fjj then Fjj{U) obeys 
a uniform U(0, 1) law and, reciprocally, if V obeys a U(0, 1) distribution then FjJ 1 ^) is F\j distributed. In 
order to use this method it is necessary to know the correlation that the random variables will have after the 
transformation, besides the function Fjj 1 . 

The method here studied consists of the following steps: 

1. propose a correlation structure for the S A held, say, the function pz A ', 

2. generate a field of independent identically distributed standard Gaussian observations; 

3. compute r, the correlation structure to be imposed to the Gaussian field from pz A , and impair it using 
the Fourier transform without altering the marginal properties; 

4. transform the correlated Gaussian field into a field of observations of identically distributed U(0, 1) random 
variables, using the cumulative distribution function of the Gaussian distribution ($); 

5. transform the uniform observations into S A outcomes, using the inverse of the cumulative distribution 
function of the $ A distribution (G _1 ). 

The function that relates pz A and r is computed using numerical tools. In principle, there are no restrictions 
on the possible roughness parameters values that can be obtained by this method, but issues related to machine 
precision must be taken into account. Another important issue is that not every desired final correlation 
structure pz A is mapped onto a feasible intermediate correlation structure r. The procedure is presented in 
detail in the next section. 
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3 Transformation Method 

Let G(-, (a, 7, n)) be the cumulative distribution function of a g^(a,7,n) distributed random variable. As 
previously stated, 

G(x, (a,7,n)) = T 2 „ _ 2q 
where T„ 1)I/2 is the cumulative distribution function of a Snedccor F distribution, i.e., 



*"■""> ~ r ft) f ft) UJ X tT l 1 + ^' '" 



The inverse of G(-, (a, 7, n)) is, therefore, 



G- 1 (t,(a )7 ,n)) = A /-lTj n 1 ,_ 2a (t) 



a 



To generate = (Z\(k, £))o<k<N-i,o<e<N-i ~ (Sa( Q! ; 1> n )i AZa) using the inversion method we define ev- 
ery coordinate of the process Z4 as a transformation of a Gaussian process ( as Z\(i,j) = G _1 (<&(C(i, j)), (a, l,n)), 
where £ = (C(*> j))o<»<JV-i,o<j<JV-i is a stochastic process such that is a standard Gaussian random vari- 

able and with correlation function (i.e. where T^((i,j), (k,£)) = E(((i,j)((k,£))) satisfying 

p ZA ((i, j), (k, £)) = Q(a , n) (T C ((i, j), (fc, £))) (3) 

for all < i, j, k, I < N — 1 and ^ (k,£) and where $ denotes the cumulative distribution function of a 
standard Gaussian random variable. 

Posed as a diagram, the method consists of the following transformations among Gaussian (INT) , Uniform (It) 
and 9^-distributed random variables: 

K II 

G- 1 

A central issue of the method is finding the correlation structure that the Gaussian field has to obey, in 
order to have the desired S A ne ld after the transformation. The function Q( a . n ) is defined on (—1, 1) by 

n (a,n)\r) \ n ) y r(n)r(-a) J 
Q(a,n)(T) = " 



l_ _ (1\ / r(n+jt)r(-a-jj) y 

fa \n) \ T(n)T(-a) ) 



1+ 



with 



R(*,n){r)= ff G- 1 ($( u ),( a ,l,n))G- 1 ($(i;),( a ,l,n)) ( /» 2 ( M ,«,T))) ( i^ 



where 

, 1 f u 2 - 2t.u.v + v 2 

<h{u,v,T) = - exp 



2^7(1— 75) "V 2(1 - T 2) 

Note that i? (a ,„)(r c ((i, j), (fc,£))) = E{Z\{i,j)Z\{k,£)) for all < i,j,k,£ < N - 1 and (i, j) ^ (fc,f)- 

The answer to the question of finding tq given pz A is equivalent to the problem of inverting the function 
Q(a,n)- This function is only available using numerical methods, an approximation that may impose restrictions 
on the use of this simulation method. 
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3.1 Inversion of Q( a ,n) 

The function Q( a , n ) nas the following properties: 

1. The set {g( a ,n)( T ) '■ T G ( — 1) 1)} is strictly included in (—1, 1), and depends on the values of a. 

2. The function Qi a>n \ is strictly increasing in (—1, 1). 

3. The values Qi a ,n){ T ) are strictly negative for all r < 0. 

Let 9( Qj „) be the inverse function of Qr a , n ). Then, in order to calculate its value for a fixed p G (—1,1), we 
have to solve the following equation in r: 

fl(a.n)T +-f-+ p-l - ^ 2 2 M =0 

v ' 1 + a \n J \ 1 (n)l (— a) J 

Then, it follows from the properties of Qf a ,n)j that for certain values of a the set of t such that this equation 
is solvable is a strict subset of (—1, 1). Table [T] shows some values of the function 9( a ,n) for specific values of p, 
n and a. Figure 0] shows r as a function of p for the n = 1 case and varying values of a, and it can be seen that 
the smaller a the closer this function is to the identity. This is sensible, since the distribution becomes more 
and more symmetric as a — > —oo and, therefore, simulating outcomes from this distribution becomes closer and 
closer to the problem of obtaining Gaussian deviates. 

Figure [5] presents the same function for a — —1.5 and varying number of looks. It is noticeable that r is far 
less sensitive to n than to a, a feature that suggests a shortcut for computing the values of Table [T] disregarding 
the dependence on n, i.e., considering r(p, a, n) ~ r(p, a, hq) for a fixed convenient uq. 

The source FORTRAN file with routines for computing the functions Qr a ,n) an d 3(a,n) can be obtained from 
the first author of this paper. 

3.2 Generation of the process £ 

The process £, that consists of spatially correlated standard Gaussian random variables, will be generated using 
a spectral technique that employs the Fourier transform. This method has computational advantages with 
respect to the direct application of a convolution filter. Again, the concern here is to define a finite process 
instead of working on 1? for the sake of simplicity. 
Consider the following sets: 



Ri 


= {{k,t) 


0<k,l< N/2}, 






= {(M) 


N/2 + 1 <k<N- 


1,0<£< N/2}, 


R s 


= {(M) 


< k < N/2, N/2 4 


1<£<N-1}, 


i?4 


= {(k,£) 


N/2 + 1 < k < N - 


1, N/2 + 1 <£<N-1}, 



R N = i?i U R 2 U R 3 U Ri = {(k,£): < k,l < N- 1}, 
R^ = {(k,£): - (N- 1) < kj< N- 1}. 

Let p: R\ — > (—1,1) be a function, extended onto i?jv by: 

(kJ)eR 2 , 
(k,£)eR 3 , 
(k,£) e Ra, 
(N-l)<k<Q<£<N-l, 
(N - 1) < £ < < k < N - 1, 
-(TV - 1) < k,£ < 0. 



p{k,£) = l 



p(N~k,£) 




if 


p(k,N-£) 




if 


p(N-k,N- 


I) 


if 


p(N + k,£) 




if 


p(k,N + £) 




if 


p(N + k,N + 




if 
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p 


a = -1.5 


a = -3.0 


a = -9.0 


n = 1 


n = 3 


n = 6 


n= 10 


n = 1 


n = 3 


n = 6 


n = 10 


n = 1 


n = 3 


n = 6 


n = 10 


-.9 




















-.953 


-.954 


-.958 


-.8 


















-.877 


-.845 


-.845 


-.848 


-.7 










-.886 


-.881 


-.901 


-.915 


-.763 


-.737 


-.737 


-.740 


-.6 










-.747 


-.745 


-.761 


-.772 


-.650 


-.630 


-.630 


-.632 


-.5 










-.613 


-.612 


-.624 


-.632 


-.539 


-.523 


-.523 


-.525 


-.4 


-.844 


-.903 


-.948 


-.972 


-.483 


-.483 


-.492 


-.498 


-.429 


-.417 


-.417 


-.419 


-.3 


-.591 


-.630 


-.656 


-.670 


-.357 


-.357 


-.363 


-.367 


-.320 


-.312 


-.312 


-.313 


-.2 


-.370 


-.392 


-.405 


-.412 


-.234 


-.235 


-.239 


-.241 


-.212 


-.207 


-.207 


-.208 


-.1 


-.174 


-.183 


-.188 


-.190 


-.116 


-.116 


-.117 


-.119 


-.105 


-.103 


-.103 


-.104 





.0 


.0 


.0 


.0 


.0 


.0 


.0 


.0 


.0 


.0 


.0 


.0 


.1 


.155 


.161 


.164 


.165 


.112 


.113 


.114 


.115 


.104 


.103 


.103 


.103 


.2 


.294 


.303 


.307 


.309 


.222 


.223 


.225 


.226 


.208 


.205 


.205 


.205 


.3 


.418 


.428 


.433 


.435 


.328 


.329 


.332 


.334 


.310 


.306 


.306 


.307 


.4 


.529 


.539 


.544 


.546 


.432 


.433 


.436 


.438 


.411 


.407 


.407 


.408 


.5 


.629 


.638 


.642 


.644 


.533 


.534 


.537 


.539 


.512 


.507 


.508 


.508 


.6 


.719 


.727 


.730 


.731 


.631 


.633 


.635 


.637 


.611 


.607 


.607 


.608 


.7 


.800 


.806 


.808 


.809 


.727 


.728 


.731 


.732 


.710 


.706 


.706 


.707 


.8 


.873 


.877 


.879 


.880 


.820 


.821 


.823 


.824 


.807 


.805 


.805 


.805 


.9 


.940 


.942 


.942 


.943 


.911 


.912 


.913 


.913 


.904 


.903 


.903 


.903 



Tabic 1: Values of function 3( Q>n ). 
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Let Za = (ZA(k,£))o<k<N-i,o<i<N-i be a S^(a, 7, n) stochastic process with correlation function pz A 
defined by 

PZ A {{kl,h),{k2j2))=p{k2-k 1 ,h-h). 

Assume that r(k,£) = 3( Q .„) (p{k, £)) is defined for all (fc,£) in i?jy. 

Let 3"(r) : i?7v — ► C be the normalized Fourier Transform of r, that is, 

N-1N-1 

y(r)(fc,£) = — ^ ^r(fc 1 ,£ 1 )exp(-2^(fc-fc 1+ £^ 1 )/7V 2 ). 

fe 1= 0^ 1= 



Let V: #jv — ► C be defined by ip{k,£) = y/^rjQ^I) and let the function : R N = {(k, I) : - (N - 1) < k,£ < 
N - 1} — > E be defined by 



0(fc, = 9^ W0(*> ^(^1^1) exp(2 7 r l (fc • h + 1 ■ h)/N 2 ), 

fc 1= £i=0 

(the normalized inverse Fourier Transform of ip) for all (k, £) € Rn', and 

( 8(N + k,£) if -(N -1) < k <0 < £ < N -1, 
6{k,£) = l 9(k,N + £) if -(N- 1) < £ < < fc < TV- 1, 
[ 9(N + k,N + £) if -(TV- 1) < k,£ < 0. 

A straightforward calculation shows that 

JV-l N-l 

{6*6)(k,£)= ^20(k u £i)0(k-ki,£-£ 1 ) = T(k,£), 
k 1= oe 1= o 

for all (k,£) e R N - 

Remark 1 We can see that 7{r){k,£) > and the last equality for all (k,£) € Rn is easily deduced from the 
results in Section 5.5 of ]12jj : more details can be seen in 1 13]. 

Finally we define C = [C(i, j))o<i<N-i,a<j<N-i by 

C(k,£) = (6*0(k,£) = NJ- 1 ((i>?(!i)))(k,£), 

where £ = (£(fc, £))(k,i).eR N is a Gaussian white noise with standard deviation 1. 

Then it is easy to prove that ( = (((i,j))o<i<N-i,o<j<N—i is a stochastic process such that Cihj) is a 
standard Gaussian random variable with correlation function tq satisfying (|3|). 

3.3 Implementation 

The results presented in previous sections were implemented using the IDL Version 5.3 Win 32 [H] development 
platform, with the following algorithm: 

Algorithm 1 Input: a < — 1, 7 > 0, n > 1 integer, p and r functions as above, then: 



1. Compute the frequency domain mask ip(k,£) — ^JJ(T)(k,£). 

2. Generate £ = {£,(k, £)) (k,e)eR N > the Gaussian white noise with zero mean and variance 1. 

3. Calculate ((k,£) = AJ- 1 ^ • ?(£)))(k, I), for every (k,£). 
I Obtain Z\{k,£) = G~ 1 ($(C(A; ! ^)), (a,l,n)), for every (k,£). 
5. Return ZA(k,£) = ^/jZ A (k,£) for every (k,£). 
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4 Simulation results 



In practice both parametric and non-parametric correlation structures are of interest. The former rely on ana- 
lytic forms for p, while the latter merely specify values for the correlation. Parametric forms for the correlation 
structure are simpler to specify, and its inference amounts to estimating a few numerical values; non-parametric 
forms do not suffer from lack of adequacy, but demand the specification (and possibly the estimation) of po- 
tentially large sets of parameters. 

In the following examples the technique presented above will be used to generate samples from both para- 
metric and non-parametric correlation structures. 

Example 1 (Parametric situation) This correlation model is very popular in applications. Consider L > 2 
an even integer, 0<a<l, <£ (for example e — 0.001 ), a < — 1 and n > 1. Let h: R — > M be defined by 



h(x) 



x if 
if 



\x\ > e, 
\x\ < e. 



Let p: R\ 



(-1, 1) be defined by p(Q, 0) = 1 if {k, i) ^ (0, 0) in Ri by 
p(k,£)-- 



h(aexp(-k 2 /L 2 )) if k>l, 
-h(aaip(-e 2 /L 2 )) if k<£. 

The image shown in Figured of size 128 x 128, was obtained assuming a — 0.4, L = 2, a 
n = 1. 



T.5, 7 = 1.0 and 



Example 2 (Mosaic) A mosaic of nine simulated fields is shown in Figure [?| Each field is of size 128 x 128 
and obeys the model presented in Example]]] with a — 0.4, 7 = 1.0, n — 1, roughness a varying in the rows 
(—1.5, —3.0 and —9.0 from top to bottom) and correlation length L varying along the columns (2, 4 and 8 from 
left to right). 

Example 3 (Non-parametric situation) The starting point is the urban area seen in Figured This 128 x 
128 pixels image is a small sample of data obtained by the E-SAR system over an urban area. The complete 
dataset was used as input for estimating the correlation structure defined by an 16 x 16 correlation matrix using 
Pearson's procedure (p below, where only values bigger than 10 -3 are shown; see appendix^). The correlation 
structure for the Gaussian process is r below, where only values bigger than 10 -3 are shown. The roughness and 
scale parameters were estimated using the moments technique. The simulated field is shown in Figure^ 



( 1.00 
0.97 
0.88 
0.76 
0.64 
0.53 
0.43 
0.36 
0.29 
0.24 
0.20 
0.16 
0.13 

\ 0.11 



0.65 
0.63 
0.58 
0.50 
0.43 
0.36 
0.30 
0.25 
0.20 
0.17 
0.13 
0.11 
0.00 
0.00 



0.22 \ 
0.22 
0.21 
0.19 
0.16 
0.14 
0.12 
0.10 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 ) 



( 1.00 
0.98 
0.93 
0.85 
0.75 
0.68 
0.56 
0.49 
0.41 
0.35 
0.29 
0.24 
0.20 
\ 0.17 



0.76 
0.74 
0.70 
0.63 
0.56 
0.49 
0.42 
0.36 
0.294 
0.25 
0.20 
0.17 
0.00 
0.00 



0.32 \ 
0.32 
0.31 
0.28 
0.24 
0.21 
0.18 
0.16 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 j 
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5 Conclusions and future work 



A method for the simulation of correlated clutter with desirable marginal law and correlation structure was 
presented. This method allows the obtainment of precise and controlled first and second order statistics, and 
can be easily implemented using standard numerical tools. 

The adequacy of the method for the simulation of several scenarios will be assessed using real data, following 
the procedure presented in Example [3J estimating the underlying correlation structure and then simulating 
fields with it. A mosaic of true and synthetic textures will be composed and made available for use in algorithm 
assessment. 
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A Estimating correlation structure with Pearson's method 

Consider the image z with M rows and N columns 

2(0,0) ••• z(N- 1,0) 

z(0,M-l) ••• z{N -1,M -1) 

and n v a positive integer smaller than min(M, N). Define n c = [N/(2n v )] and nt = [M/ (2n v )], where [x] = 
max {k 6 N : k < x} for every real number x. For each i = 0, . . . , n c — 1 and each j = 0, . . . , rif — 1 define c(i,j) 
the submatrix of z of size 2n v x 2n v given by 



z(2n v i, 2n v j) 



z(2n v i + 2n v — 1, 2n v j) 



z(2n v i, 2n v j + 2n v — 1) • • • z(2n v i + 2n v — 1, 2n v j + 2n v — 1) 
and let z v (i,j) be the submatrix of c(i,j) of size n v x n v given by 



zv(i,j) 



z(2n v i, 2n v j) 



z(2n v i + n v — 1, 2n v j) 



z(2n v i, 2n v j + n v — 1) • • • z(2n v i + n v — 1, 2n v j + n v — 1) 



We will consider that z v (i,j), for every i = 0, . . . ,n c — 1 and every j — 0, . . . , rif — 1 is a sample of the 
random matrix 



Z(0,0) 



Z(n v - 1,0) 



Z(0, n v — 1) ■ ■ ■ Z(n v — 1, n v — 1) 
The autocorrelation function of the random matrix Z is defined as 

E(Z(m, n)Z(k, £)) - fj,z(m, n)fi z (k, £) 



Pz((m,n), (k,£)) 



az(m, n)a z {k,£) 



where fi z (k,£) = E(Z(k,£)) and a z (k,£) = y/Var(Z(k,£)), for every 0<m,n,k,£< n v - 1. 

The function pz can be estimated using Pearson's sample correlation coefficient based on z v (i,j), i 
0, . . . , n c — 1 and j = 0, . . . , n/ — 1, i.e., for < m, n, k, £ < n v — 1 by 



r z ((m,n), (k,£)) 



C z ((m,n),(fc,£)) 
sz(m,n)s z (k,£) ' 
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where 



n.f-1 n c -l 

Cz((m, n), (k,£)) = ^ ^ (z(2n v i + to, 2n v j + n) — z(m, n)) (z{2n v i + fc, 2n v j + i) — ~z(k, £)) , 

j=0 i=0 



nf—1 n c —l 



, ^ (z(2n v i + to, 2n„j + n) - z(m, n)Y 



s z (to, rc) 

z(m, n) = ^ ^ z(2n v i + to, 2n v j + n). 



Uf-l n c -l 



n c nt 

J j=0 i=o 
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Figure 6: Correlated S°(— 1.5, 1, redistributed amplitude image with the correlation structure defined in Ex- 
ample [TJ 
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Figure 7: Mosaic of nine simulated fields. 
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Figure 8: Urban area as seen by the E-SAR system. 
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Figure 9: Simulated urban area using a non-parametric correlation structure. 
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